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The original sandpile model of Bak, Tang and Wiesenfeld from 1987 has in- 
spired lots of consequent work and further ideas of how to describe the birth of 
scale-invariant statistics in various systems and in particular models. In this article 
the basic ingredients of self-organized criticality (SOC) are overviewed. In line 
with the orginal arguments of BTW SOC is now known to be a property of systems 
where dissipation and external drive maintain a delicate balance. The qualitative 
and quantitive understanding of the SOC state and deviations from it can thus be 
addressed, by mapping the typical cellular automata exhibiting SOC to theories, 
that in fact describe critical extended systems as such. Currently two such ap- 
proaches are known, based on the connections of SOC to variants of absorbing 
state phase transitions and the physics of elastic manifolds or interfaces in random 
media. These are reviewed and discussed. Finally, some open theoretical problems 
and experimental suggestions are outlined. 

1 Introduction 

Usually statistical physics is on the textbook level a subject that deals with "bor- 
ing" properties of matter and the associated basic thermodynamic quantities. There 
is however a modern challenge to this, that attempts to deal with all kinds of phe- 
nomena with statistical physics tools. One of the simple reasons is the discovery 
of many (effective) power-laws or "fractals", in the statistics of measurable quanti- 
ties. This implies that within the appropriate window the system is scale-invariant, 
which in the language of statistical mechanics brings us to the field of phase transi- 
tions. The next question that follows is then: how come so many phenomena have 
parameters finely tuned to the point where the equations and laws governing them 
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allow for "criticality"? Consider the prototypical two-dimensional Ising model as a 
example: for sure one needs to tune the temperature to Tc to reach the same effect. 

The seminal papers of Bak, Tang and Wiesenfeld (BTW) started a literal avalanche 
in this respect. They coined the term of "Self-Organized Criticality" (SOC), with 
an associated set of claims, conjectures, results, and a simple model well within 
grasp of anyone fluent with computers. The main idea of BTW was that such a 
sandpile mode would be effective in explaining the presence of 1// -noise in Na- 
ture, a prototypical example of power-laws indeed. The purpose of this article is 
to give an overview of where the mainstream SOC understanding has lead to, from 
the starting point of the 1987-1988 BTW publications [1|. For this goal it is in- 
structive to start with the original BTW sandpile model, for a brief look as to what 
the rules and aspects of the SOC sandpiles imply. These are often most defined as 
cellular automata, on hypercubic lattices of size . The original BTW numerical 
results concentrated on d = 2. Each site x has z{x, t) grains. There are two crucial 
ingredients that define the dynamics in the BTW model. When z{x,t) exceeds a 
critical threshold Zc (a constant, here 3), the site is active and topples. The grains 
(4 in the 2d BTW example) are removed from x and given to the nearest neighbors 
(nn). The effect of this rule is to create a non-linearity: in the absence of any activ- 
ity in the nn's of x the site stays quiescent, forever if Zx = z{x, t) < ■ In other 
words, the BTW model has Zc equal to a constant, Zc = 2d — 1, and the toppling 
rule: 



where y denotes all the 2d nearest neighbors of the site x. 

One may now study two fundamental limits (note that clearly other scenarios 
can be envisioned). Either one prepares the system with a fixed number of grains 
ntot = Yl,x , and observes what happens, with e.g. periodic boundary conditions. 
Due to the presence of the nonlinearity it is easy to see that there are two opposite 
limits, an absorbing state where all Zx < z^ and nothing happens, and a state of 
eternal activity where for all t for some x Zx > Zc ■ This defines a "fixed energy 
ensemble" (FES) as it has been denoted in the literature. Clearly there has to be a 
(phase) transition as n tot is varied between the two extremes. 

However, the particular case of "original" SOC of BTW is created by different 
conditions. What one wants to obtain is a steady state. This is obtained by com- 
bining open boundary conditions that are balanced with a "drive". The boundaries 
are chosen to be such that the grains which "topple" out of the system are lost, 
simply. Eg. in 2d a site which touches the boundary loses a grain out of the four it 
gives out. The SOC state is now obtained by using a modern version of Maxwell's 



z{x,t + l) = 
z{y,t + l) 



z{x,t) - 2d 
= z{y,t) + l 



(1) 
(2) 
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demon: if and only if tiiere are no active sites {z^ > Zc) one grain is added to a 
randomly chosen site x: z{x, t) — > z{x, t) + 1. Except for the presence thereof 
there is no "tuning" whatsoever. Clearly, the effect of an added grain is either to 
build up Zx or to launch an avalanche if x is "marginally stable" if Zx = Zc ■ The 
original claim of BTW was that the sizes of avalanches - measured as the number 
of times that any single site would get active after a single grain addition - would 
follow a power-law, as would the duration (multiple simultaneous topplings being 
allows) and the area or support on the L"^ lattice, 

f{s) ^ s--^f{s/s,), (3) 

with the avalanche size exponent Tg and the cut-off dimension Sc = L^'' ■ The 
scaling function / expresses the fact that avalanches are limited by the scale Sc and 
thus decays quickly beyond that. Thus there would be a subtle balance of driving 
and dissipation in a SOC system, with power-laws arising apparently without any 
fine-tuning and hence the birth of the acronym. 

The essence of SOC as from the BTW model is thus non-linearity and a combi- 
nation of dissipation and driving |2|. There is an easy set of questions that follows 
immediately, in particular from the analogies with other phase transitions: 

• universality: what kind of exponents (r^ , Ds etc.) would one be able to 
obtain depending on the particular details of the model? 

• what does SOC "mean" in fact? 

• how broad is this paradigm in terms of applicability? 

• what is the right "continuum theory" for SOC? 

In the rest of this article the status of the two ideas and the follow-up ques- 
tions will be examined, in the light of recent progress on mapping SOC sandpiles 
to other systems. In this, the main emphasis is on analogies with interface de- 
pinning and absorbing state phase transitions. Both these allow to adapt to the 
ensembles (SOC, FES) in which the simple SOC automata are run, and moreover 
present the essential feature of a "nonlinearity", of a threshold element as we will 
see below. These present two different "generic" types of continuum limits. For 
the second one, we have as the dynamical variable the activity, p(x, t). The in- 
terface depinning variants of SOC boil down in most effective form as a variable 
ff(x, t) = /* /3(x, t)dT and its dynamical equation dtH{x, t) . . .. The picture of 
SOC is completely equivalent as to how the right ensemble is created. This will be 
highlighted with examples in the next section. 



3 



1.1 Theoretical aspecs of sandpiles 

Analogies with other systems in statistical mechanics have been around for a long 
while but have not been exhausted even by now. Tang and Bak already noted that 
one can perturb - on the cellular automaton level - sandpiles by driving them at a 
steady, slow rate so that overlapping avalanches are not too frequent. They also 
conjectured - based on ideas from equilibrium phase transitions - several exponent 
relations, to describe the correlation length and the correlation time in the prox- 
imity of the vanishing drive rate, and most importantly identified a (correct) order 
parameter, the average activity or probability that a site is active {zx > Zc in the 
BTW case). The mean-field treatment results in the exponents of the contact pro- 
cess or directed percolation (DP) |3l, and hence it immediately brings to the ques- 
tion whether sandpiles are just an example of an absorbing state phase transition 
like DP. 

For the BTW model, again, Narayan and Middleton noted that it can be ob- 
tained from a discrete version of an equation describing charge-density waves (not 
surprising considering the history of the model as such) |4l. This equation is 
known another contexts and with a different noise term as the quenched Edwards- 
Wilkinson equation (QEW). It describes the dynamics of a forced domain wall in 
a random magnet. The impurities pin the interface and a driving force counteracts 
the pinning. Later, Paczuski and Bottcher pointed out that a particular "rice pile" 
SOC model |5 1 could be mapped also into the same system. 

The presence of these two analogies (see also [6.2.1) or in more concrete terms 
mappings to (continuum) systems from SOC ||8| |9| makes it possible to identify 
firmly many answers to the above mentioned problems, and in general to under- 
stand the critical properties without any conceptual confusions. These are dis- 
cussed, in combination with the missing pieces of the puzzle, in the later para- 
graphs. As will become clear the nature of SOC in such typical cellular automata 
is a particular ensemble version of various non-equilibrium phase transitions. In 
many of the early papers analogies were drawn with the physics of equilibrium 
ones. Here, however, the situation is different and one for instance is not nec- 
essarily able to define a free energy for the system, or even a Hamiltonian (e.g. 
the perennial Kardar-Parisi-Zhang problem serves as an example OSl)- It also be- 
comes more evident what the correct renormalization group procedure is, in con- 
trast to various real-space, mean-field and other attempts 111, 12. 13. 141. This 
leads to an understanding of where the differences in the numerically found scal- 
ing behavior among the BTW, Manna 1(151. Zhang [ 16 1 etc. models originate from. 

In most cases it is also obvious that the SOC state is indeed achieved only 
in a particular point of the phase space, the only reason why this is not obvious 
to begin with is that it becomes apparent after such connections to other models 
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are made clear. Thus in general there is no "generic scale invariance" lITTl ITSi 
unless one goes one step further: in the driven interface paradigm the system still 
exhibits critical fluctuations off the SOC state, but the character of the fluctuations 
changes. The simple analogy here is the thermalization of the noise in a depinning 
problem, and the changed interface fluctuations. Here one has to be careful about 
the ensemble that is perturbed off the critical point: in translation-invariant models 
with a uniform dissipation rate e one can observe avalanches when e — > 0^ fT3llT4l 
similarly to the SOC ensemble. This is of course true only when the drive rate h 
is adjusted so that h/e fr7lfT9lfT4ll . ie. in the limit that these two very slow 
scales are still "infinitely separate" (Tt\ . If one accepts the idea that SOC actually 
arises from "normal criticality" by some clever mechanism (e.g. "sweeping the 
instability" as was proposed |20|), then the question still remains how this can be 
made to work, and in e.g. the interface picture this becomes immediately obvious. 

The original BTW model itself still attracts interest since it does have a power- 
ful Abelian symmetry, peculiar of a discrete cellular automaton. As shown by Dhar, 
Priezzhev and co-workers in a series of papers there are connections to spanning 
trees on lattices [21 , 22 ], thus to the critical q = I Potts model (and in 2d to confor- 
mal invariance) 1231 . In this respect it has however lost the aura of generality that 
the original papers advocated, not forgetting the long-standing controversies about 
the various universality classes of SOC models - related to the numerically deter- 
mined exponents based on Eq. (3). Recent numerical studies have even pointed out 
that the BTW model shows multiscaUng: the avalanches exhibit a full spectrum of 
scales, whose origins are not understood |24|. Such numerical studies have made it 
apparent that there are definately different universality classes, as the connections 
to interface depinning and APT's should make clear below. 

In the following models that involve strict "extremal dynamics", in the sprit 
of the Bak-Sneppen model |25|, are excluded. The idea here is that one (in the 
interface context) advances the site which has the largest local force. Clearly this 
kind of dynamical mechanism is just another version of the demon inside the BTW 
sandpile, an even smarter one since he has now to know the state of the system in 
greater detail. We shall also not discuss the so-called forest-fire models ['26', IIT7l . 
that have the difference that one has three states for a site: empty, a living tree, and 
a burning one. For this class one has no such conservation law as exemplified by 
the Laplacian of the interface equation or by the energy field of the APT class with 
conservation. It is an open question how to find similar connections to continuum 
models. We also omit directed models which as such are often more simple to 
understand theoretically. 
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1.2 Experimental signatures of SOC 

The obvious question to ask is, what is the quantity to measure given that SOC is 
in fact a "usual phase transition"? There are two answers: either one is concerned 
with avalanche behavior, hke in a typical sandpile model, or then one is interested 
in some quantity that would be able e.g. to establish whether the dynamics are 
indeed governed by or close to something akin a SOC critical point. Given the 
fact that one needs a separation of timescales it becomes clear that sandpile-like 
behavior is not easy to conjure in an experiment. On the other hand, it is clear that 
a generic phenomenon like l//-noise is precisely speaking far from the typical 
SOC manifestations. 

Laboratory-scale experiments for this purpose were attempted soon after the 
first BTW papers inspired by the "sandpile" idea l28ll . Unfortunately, there is 
little reason as to why a real sandpile should behave like a cellular automaton. In 
particular, the flow of grains once started involves inertia 1.29 J . The ingenuity of 
experimentalists then lead (in Oslo) to the discovery that one can simply substitute 
with (elongated) rice grains, and a power-law avalanche distribution foUows l30l 
(see also |31 1). 

In essence, to mimick sandpiles one thus needs diffusive motion of the "par- 
ticles", ie. the equation of motion has to contain precisely the same ingredients 
of a threshold and overdamped motion once that is exceeded. Such conditions are 
provided by domain walls (DW) in ferromagnets - as we will see one can easily 
define a sandpile that describes the dynamics of a DW - and vortices in type-II 
superconductors. The former can be observed by e.g. optical means or by the 
noise produced by domain wall motion, known as the Barkhausen effect l32ll . It 
is possible to observe avalanches, with a size distribution that extends as a power- 
law over more than three decades I33ll34. ,35„ 36, .37,1 . The classical Barkhausen 
DW experiment corresponds however to the slowly driven, translationally invariant 
case outlined above (and not to the SOC ensemble): a constant external field ramp 
corresponds to a constant velocity for the driven DW. 

The recipe for superconductors is to increase the external field which pushes 
vortices or flux lines into a sample (recall of type II). The Bean state ll38l l391 that 
ensues has a conservation law (vortices), and can be controlled with the rate of 
increase of the external field. One can observe avalanches, exhibiting power-law 
behavior 1.40. . 4 L . 42 1. Indeed there are sandpile models to describe such behavior 
ll43l . and on the other hand coarse-grained lattice models that are close to sand- 
piles in spirit l44l (to say nothing about more complicated simulations [45 1). In 
both these cases the actual role of SOC in the phenomena at hand seems slightly 
secondary, however. 

Avalanches have also been observed in a more promising field from the point 
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of view of potential applications: fracture and plasticity. While there is lots of 
evidence (in terms of acoustic emission, the energy relased by microscopic events) 
for the existence of power-laws both in laboratory experiments f46'. El and in 
earthquakes, it is not clear at all how much these have to do with SOC. In failure 
of brittle materials it is clear that the underlying processes are irreversible and have 
nothing in common with the basic premises of a SOC state. In the case of earth- 
quakes the scales involved are such that it may seem more prudent to explain eg. 
the classical Richter's and Omori's laws (of magnitudes and aftershock intervals) 
via a suitable SOC model. However, no direct connection has been established so 
far B9l . Plasticity and dislocation dynamics seem better in this respect, as a field 
of application of SOC ideas (see l50ll5T]| for AE evidence). One has the possiblity 
of a steady-state and the system can be driven by slowly applied external strain. 

Perhaps the most enticing field where SOC concepts have been utilized is 
plasma and space physics. The crucial ingredients here are the presence of non- 
linearities (given by the character of e.g. the magnetohydrodynamic equations that 
can be proposed for the Earth's magnetotail sheet), external drive (ionospheric ac- 
tivity vs. solar wind, or the solar flares vs. sun's heating), and spatial inhomogene- 
ity in that the dissipation and drive do not necessarily take place uniformly. Both in 
space and laboratory plasmas one can obtain measured statistics (like the particle 
flux driven by the nonlinear transport at the plasma edge in a fusion device, or the 
magnetospheric AE index) that can be then analyzed in terms of SOC character- 
istics. Several reviews exist that outline the subject (whether one is interested in 
fusion or space aspects) I52ll53ll54ll55ll56ll57ll . There is clear evidence of power- 
laws in eg. the waiting times of avalanches and energy dissipated |53|. In these 
contexts some of the differences to usual sandpile models are based e.g. on com- 
parisons to (thresholded) real signals, and on the presence of correlated and varying 
drive signals. Both these make comparisons between theory and experimental or 
observational quantities somewhat hard. 

1.3 Overview 

Next we start by presenting a pedestrian picture of the connection of SOC sandpile 
models to those of non-equilibrium statistical mechanics. The following section 
goes on to a discussion of such mappings in more detail: the universality class 
of absorbing state phase transitions with a conserved field is outlined as the first 
one. Section IV considers SOC in the light of interface depinning. First the basic 
background about the latter is given. Then the mapping of sandpiles to discrete 
interface equations is discussed, together with the implications that the various 
terms and the ensemble(s) have. To complete the picture we also consider the 
possibhty of extending the known SOC ensembles, by the quenched KPZ equation. 
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Finally, in the conclusions we List a number of open topics for further research, and 
summarize the state of the art. 

2 Two models and basic ideas 

To illuminate the differences that arise from sandpile rules and the connections to 
non-equilibrium phase transitions we next discuss some examples. Consider for 
this sake two models, both defined in one dimension, on a lattice with x = I . . . L. 
In the first one, the rules are similar to the BTW model. Each site has grains. 
For simplicity, we start with the FES ensemble and prepare the system with one 
grain per site. The difference to the BTW one is that now the thresholds Zc{x) are 
chosen to be random. We assign to each x a. Zc value of 2 with a probablity 1 — p 
and with probability p after each toppling. In this way, Zc{x) is an iid random 
variable, if one considers its values at two x, xq or at the same site, separated by 
toppling(s) at x. Then the toppling rule is for Zc < z^ 

z{x,t + l) = z{x,t) -2d (4) 
z{y,t + l) = z{y,t) + l (5) 

where y are the nn-sites of x. The dynamics of the model is simple: one compares 
the number of grains at a site to the local threshold. In the example above, the 
Zx follows simply from differences in the flux from neighbors (one grain per nn- 
toppling) and the flux out due to local topplings, One can define a 'local force' 
as 

f{x, t) = riin - Uout - Zc{x) (6) 

in terms of nj„ (grains added to site x up to time t) and Uout (grains removed from 
x). The n-variables can be directly be interpreted in terms of an interface variable, 
H{x,t), or history one which follows the memory of all the activity at x. 

In particular, the above sandpile turns out to be the Leschhom-Tang (LT) cellu- 
lar automaton 1 58 1, used to simulate interface depinning for the Quenched Edwards- 
Wilkinson (QEW, or Lineal- Interface Model, LIM) equation ||59l|60l. The LT au- 
tomaton follows an interface H{x) at each discrete time step tj with the equation: 

H{x, ti+i) = H{x, ti) + 1, /(x, ti) > (7) 
= H{x,ti)J{x,ti)<Q (8) 

where the force / is in the QEW/LIM language a combination of elasticity and a 
random quenched pinning force 

/(x,t,) = V2F(x,t,) + ??(x,i?). (9) 
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H (change) 
H (before) 
H (much later) 



avalanche, dH(x) 



X 



H(0) = 



H(L+1) = 



Figure 1 : One-dimensional schematic example of the interface or history represen- 
tation of a SOC model. The mean-field interface is parabolic, which implies that in 
the SOC steady-state v{x) is parabolic, too. Notice the boundary conditions H = 
that ensure the loss of particles (equalling the increased elastic energy). 

V^if(x) = (Id) i^j+i — 2Hi + Hi^i is the discrete Laplacian. The noise in the 
original LIM cellular automaton is 



This choice implies that there is a average driving force F = (/) = 2p — 1, which 
is the control parameter (the brackets denote ensemble averaging). The critical 
point is estimated to be at pc ^ 0.800 |58|. Note in particular that (obvious) fact 
that the FES critical point is that of the QEW universality class, to which we return 
later. 

So we have discovered that we can easily associate a sandpile model with the 
QEW cellular automaton. If one would use open boundary conditions this would 
mean that in the SOC ensemble, with no active sites, one grain is added to a ran- 
domly chosen site, z{x, t) z{x, t) + 1. One can now look at f{x, t) in this case, 
and notice that there is a "columnar" force term F{x, t) which counts the number 
of grains added to site x by the external drive. In the original example F = 1. In 
this ensemble the non-linearity is provided by the random thresholds inside /, and 
the slow drive is the same as in the QEW case, basically. 

The SOC ensemble works (Figure (1)) so that the local force is increased in 
steps, by adding a grain at x and making F{x,t) — > F{x,t) + 1. Sometimes, 
/ > and an avalanche starts. The interface moves at x by one step, 5H{x) = 1. 
In the subsequent dynamics during the avalanche the columnar force term F{x) 
does not change. The right choice of the interface boundary condition is = 
which is to be imposed at two "extra sites" (x = 0, x = L + 1 for a system of size 




(10) 
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L in Id). The {F{t)) increases slowly with time, so that the shape of the interface 
is on average a parabola. There are now three easy-to-define questions following 
from the SOC setup: i) can one understand the "avalanches" (bursts of activity) 
from the QEW depinning, with its presumably well-known exponents? ii) what is 
the effect of the parabolic interface shape on the physics of depinning, if any? (this 
relates to i) if things get complicated) and iii) how does the dynamics of such a 
driven parabola compare again to that of the "usual" QEW, if it is perturbed off the 
critical point? 

The other model is simply constructed by considering activated random walk- 
ers. One puts Zx of these per site, and the rules are such that if z^: = 1 the waUcer is 
frozen, and for Zx > 1 per unit time two walkers take off from x. This is in actual 
sandpile terms the Manna model 1 15 1, if the walkers are called "grains". We have 
that Zc = 1, and the toppling rule 



with y' and y" two randomly chosen neighbor sites of x. 

In the FES case (consider periodic boundary conditions) there is naturally a 
phase transition located at ntot < L. Like in the QEW, the bare value of the con- 
trol parameter (here ntot / L) at the critical point is hard to compute. The important 
point is that the model is clearly an example of an "absorbing state" phase tran- 
sition: any product measure of states with z^ = 1 or will do. Thus there is an 
infinite number of absorbing states, and moreover trivially ntot is conserved. These 
suffice to make it possible - and clear - that there can be differences to e.g. the con- 
tact process (this would amount in the random walker language to a model where 
walkers can die out, and have off-spring, moreover Zc = 0), in the critical behavior. 

In this model there is as well a memory effect: any configuration of grains in a 
local neighborhood is frozen till an avalanche sweeps over it readjusting the Zx's. 
Note that the effect of this is diffusive since if a grain moves around, it does so by a 
simple diffusion process. The SOC ensemble is again obtained similarly to the LT 
automaton by allowing walkers to escape at the boundaries, and adding a walker at 
a time randomly (zx —>■ z^ + I) if no active sites exist. In this case, it is easy (and 
easier than in the QEW case) see why there might be subtle differences between 
the two ensembles. The average {zx) in the quiescent state between avalanches 
is created by the average flux of grains out of the system. In ID, clearly any site 
for which x ^ L/2 there has to be a net flux of grains towards the absorbing 
boundaries. Thus the profile of (zx) will for sure have a gradient as a function of 
X, giving rise to "finite size" effects that are difficult to analyze. 



z{x,t + l) 
z{y',t + l) 



z{x,t)-2, 
z{y',t) + l, 
z{y",t) + l 



(11) 
(12) 
(13) 
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3 Sandpiles and absorbing states 



The classical model for a system with an absorbing state is the "contact process", 
where particles diffuse, die and are born as off-spring from neighboring ones on a 
lattice. In terms of a coarse-grained density-field p one can illustrate its behavior 
with the mean-field equation 

f = (A - 1)/. - (14) 

that includes the competition between the two mechanisms affecting the density p 
ll3l l6Tl . The MF- variant clearly has a phase transition at a Ac = 1, with p = 1 — 1/ \ 
in the stationary state. The physics version on a lattice is directed percolation, 
which has as well an upper critical dimension du = 4. The DP can be analyzed 
field-theoretically with standard coarse-graining techniques, and there is an exten- 
sive body of work on the critical exponents and scaling functions etc. for d < ci„. 

The exponents that are of interest here are such that they describe the behaviorat 
the critical point and in its proximity. The former gives rise to u±, i/y, z, and /?. 
These describe in turn the temporal and spatial correlation scales, the dynamics of 
the correlations, and the order parameter for A > Ac. One has the scaling relation 

p{A,L) = L-'^/''^R{L'^/''^), (15) 

with A = A — Ac the distance to the critical point, i? is a scaling function with 
R{x) ~ for large x. For L » ^ ~ A^'^^ we expect p ~ using for 
the correlation length. When A = we have that p{0,L) ~ L^f^/'^^. Above the 
critical point the order parameter has a stationary value. 

Such exponents and the scaling contained in Eq. (fT5l would then be the goal for 
an absorbing-state phase transition description of sandpiles. The example of acti- 
vated random walkers aka the Manna model of the previous section makes it rather 
obvious that the FES ensemble - with the mild caveat that the transition should be 
continuous - might be described analogously to the DP. The Manna universality 
class is not expected to be in the DP universality class, due to the presence of the 
conservation law and the presence of an infinite number of absorbing states. The 
questions that remain are then highly non-trivial: first, what is the right effective 
field theory? Second, are there complications related to the specific character of 
the SOC ensemble? 

The former of these two questions can be answered exactly for a class of 
stochastic models, with an infinite number of absorbing states, and with a static 
conserved field (NDCF, non-diffusive conserved field class). It is the coupling of 
this field to the order parameter evolution that creates a unique universality class 
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Steady state exponents d = 2 





P 




f3/iy± 


z 


e 


CTTP 


0.64(1) 


0.82(3) 


0.78(3) 


1.55(5) 


0.43(1) 


CRD 


0.65(1) 


0.83(3) 


0.78(2) 


1.55(5) 


0.49(1) 


Manna 


0.64(1) 


0.82(3) 


0.78(2) 


1.57(4) 


0.42(1) 


DP 


0.583(4) 


0.733(4) 


0.80(1) 


1.766(2) 


0.451(1) 



Table 1: Critical exponents for steady state simulations in d = 2. Models: CTTP: 
Conserved threshold transfer process, CRD: Conserved reaction-diffusion model, 
Manna: Manna sandpile, DP: Directed percolation. From 1641 . 



11621 1631 1641 . This is a conjecture supported by extensive simulations of a number 
different models ||62l|63l: a conserved threshold transfer process l65l . a conserved 
lattice gas with repulsion of nearest neighbor particles, and a reaction-diffusion 
model (CRD) with two species of particles, A and B. The exponents (in 2d ) are 
illustrated in Table I, together with the CP ones. 

The conserved reaction-diffusion model studied in ll63l has the pleasant fea- 
ture that it can be mapped exactly, using a Fock-space representation and creation- 
annihilation operators [66 67 1, into an effective action, or equivalently into a set of 
Langevin equations |63|. The resulting theory has, forgetting some naively irrel- 
evant terms upon power-counting — , a structure proposed earlier on phenomeno- 
logical grounds ll68l . This set of equations reads 



dt 



= ap{:i^) -bp{ji.y + V^p{jc,t) - nilj{:!c,t)p{:)c,t) (16) 
+ cr\/pixv,t))r]{x,t) 



DV^pi^,t) (17) 

For p, this looks like a Reggeon field theory 1 69 1 (used to describe DP) which is 
coupled to an extra conserved non-diffusive field, t/j. One thing is immediately 
clear, that is the theory contains a term that accounts for the memory of the dynam- 
ics of the grains in the model (0, again) which then couples to the actual activity 
p. These fields ai^e also often called "energy" (grains) and "activity" (active sites), 
respectively. Due to the presence of this coupling the theory is non-Markovian. 

Above, ry is a Gaussian white noise with a trivial correlator (?7(x, t)r/(x', t')) ~ 
(5(x — x')5(t — t'). The parameters (a and so forth) are fixed. Notice the presence 
of the linear term, that orignates from the initial configuration. The effective noise 
strength is linear in the local density, since in the coarse-graining the activation and 
passivation of particles on the microscopic level becomes a Poissonian variable, 
with the variance equalling the expectation. In particular it vanishes for p = 0. 
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The second equation describes how the conserved density diffuses upon the 
influence of actual activity. In the absence of such the dynamics is frozen. It can 
be integrated out formally, to obtain a single equation for the activity p. More 
concretely 



The first contribution in Eq. ( fT^ follows from the initial condition, and is a 
"quenched" term: it is frozen and not affected by the presence of the thermal noise. 
The second is a non-Markovian term. The theory as it stands has one easy conse- 
quence: du = ^ (though there is a claim by Wijland that in fact = 6 iT/Oll ') in 
agreement with numerics on this universality class. The problem is that Equation 
(IT6l is difficult to renormalize due to the presence of the extra terms (in addition to 
the DP-RET) ll64l . Thus one has no useful predictions as such. In the next section 
we discuss the relation of SOC to depinning, in which case the corresponding sit- 
uation in the QEW becomes pertinent. Eor quenched depinning the lesson is that 
the quenched noise field renormalizes in a highly non-trivial fashion (in particular 
from the viewpoint of the technicalities of the computations). Here an analogy 
would be the correlated activity that reflects the landscape (grains or particles), at 
various instances of time. 

It should also be noted that the theory of Eq. (fT6l fails in the presence of a 
number of effects - or does not hold of course for all possible sandpile models. The 
BTW model has no randomness in the toppling rules, its conservation laws play 
a role and the situation is different |21 1. Eor the LIM automaton of the previous 
section, the additional quenched noise clearly plays a role in the way the activity 
stops (or is "pinned"). The question what the NDCE class means as a depinning 
problem will be discussed in the next section, but it is worth underlining that the 
connection between the two kinds of transitions would merit further study |71 1. 

The SOC state in the EES case arises on a simplistic level by combining an 
absorbing state with another demon (a close cousin of the BTW one), with the task 
of driving the system if p = 0, everywhere. This is combined by a dissipation (like 
the boundary losses) that ensures that for p > dp/dt < 0. This is of course just a 
complicated way of stating the simple rules of CA's (since we know that the state 
is such that {p) = 0"'" due to the separation of timescales with respect to dissipation 
and driving). Technically this implies that (z^) in the SOC case is limited above 
by the critical value of the EES critical point. This is trivially so since the opposite 
would imply that {px) > 0, locally for some x. The actual scaling function of 
{zx){'K, L) is however not known. This is in fact one of the fundamental questions 
of SOC sandpiles: if one can separate "bulk" behavior from boundary effects. 

Recent work by Dickman has shown that precisely at the SOC critical point 
the sandpile avalanche properties obtain logarithmic corrections to their scaling 




(18) 
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functions with L. This surprising finding implies that the finite size scaling in 
the SOC ensemble is not easy to understand based on the translationally invariant 
FES (or depinning case). This would also be expected to hold in general for the 
NDCF-class (for all its models in the SOC ensemble). 

Notice however again that the description in terms of the coarse-grained the- 
ory invites for two other scenarios: clearly the critical point can be approached by 
any suitable dissipation mechanism, not only the boundary loss -one. A possibil- 
ity is bulk dissipation (e.g. a constant probability per step that a diffusing grain 
is removed). Such perturbations are closer to the FES case, since translational 
invariance is restored. 

4 SOC models and depinning 
4.1 Interfaces in random media 

Rough interfaces fT^ fTnUT^rT^ are objects with "self-affine" geometry. They can 
move due to an applied force, such as a magnetic field or a pressure difference. The 
essential point is that the thermal fluctuations may be neglected, if the noise is due 
to the presence of impurities or defects in the material. Such disorder is quenched, 
static in time. The interface can get locally pinned, if a local "pin" is strong enough 
to overcome (generalize) surface tension. A characteristic example is the quenched 
Edwards-Wilkinson (QEW), or linear interface model (LIM) 1591 1601 l58ll : 

a IT 

— = S/^H + r]{x,H)+F. (19) 

Here H is driven by F, and the combination of the surface tension and the ran- 
domness r] gives rise to the critical behavior In particular, depinning takes place 
at a force F close to a critical value Fc. Depinned interfaces move with a velocity 
(order parameter) u ~ with f = F — Fc > 0. Pinned interfaces are blocked by 
"pinning paths/manifolds" which arise from the quenched disorder in the environ- 
ment. 

The description of kinetic roughening involves analogously to the absorbing 
state phase transitions a correlation length ^. Statistical scale invariance assumes 
(statistical) translational invariance, developing both in time and space f72l ITOl . 
The typical quantity to measure is the roughness (mean square fluctuation) Wq{L, t), 
where L is the measurement scale. One often observes power law scaling for its 
moments, 

Wq{L, t) = (((5i?(x)5))i/^ _ ^x, for L < ^{t) (20) 

saturating to a constant for L larger than ^ (the inner brackets imply averaging over 
X, and the outer an ensemble average). Xg is called q-th order roughness exponent. 
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: Interface moves over an area 



Figure 2: An interface in a random enviroment: a part moves under the influence 
of F, the rest is pinned. 



allowing for multiscaling should the Xq differ from each other (as for "turbulent 
interfaces" [ 75 1, when the correlations are dominated by the largest height differ- 
ence between two neighboring points, /S.H = \Hx+i — Hx\ 17511761 177I '). 

Usual self-affine scaling means Xq = X = const, and ^{t) increases as ^{t) ~ 
t^/^, with z the dynamical exponent. The maximal extent of interface fluctuations 

w{t) = lim w{L, t) ~ i{tY ~ t>^'' = t^, (21) 

L— >oo 

which defines the exponent /?, combining into a scaling form 

w{L,t)=a{t)atr^l^-^)- (22) 

with the scaling function ~ for x < 1 and a constant for x » 1. If the 
amplitude a{t) increases, so-called anomalous scaling ensues EH, and if X > 1 
"superroughness" |78l|79l. The LIM in Id is an example ll58l ISOl ISTI 821. and 
below we will discover that several SOC automata in the interface description have 
this property as well, x can also be described by the structure factor, or spatial 
power spectrum 

5(k,t) = (|/.(k,t)| ) ^ I ^ ^ ^^^^^^ (23) 

where /i(k, t) denotes the spatial (d-dimensional) Fourier transform of /i(x, t). 
Also, the height difference correlation function G2(x, t) = (|/i(x, t)— /i(0, = 
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^(k, t) {1 — cos(k • x)) can be used, but note that it reflects only the local rough- 
ness exponent x\oc and thus for superrough interface, with x > 1. one needs S (for 
L fixed). For temporal quantities the analogue is naturally the power spectrum, in 
the stationary state, S{'k,uj) = (|/i(x, (which has obvious problems in the 
SOC state, due to the separation of timescales), or correlation function of the gth 
moment of height "jumps" over a temporal distance t 

Cg{^,t) = (|/i(x,t + s)-Mx,s)r)i/'?. (24) 

Generally one finds an increase Cq(x, t) ~ over "short" time distances t. In the 
case of (conventional) scaling it is related to the early time increase of the width, 
g-moments Pq = (3 = x/ z. It is easier to use than 'w{t) at early times. 

The main universality class of the QEW, Eq. (fT9l . is reached for random bond 
and random field disorder (these flow in the RG into the same fixed point). The RF 
noise correlator reads 

{r]{x, h{x, t))T]{x', h{x', t'))) = A{h{x, t) - h{x' , t'))6{x - x') (25) 

where A{u) is in practice taken as a delta function, A{h{x,t) — h{x',t')) 
6{h{x, t) — h{x', t')). If the velocity v is finite, the average motion and the fluctu- 
ations separates, hence for Ifh > {vty^^, 

6(h(x, t) - h(x', t')) 6(v(t - t') + 5h) -5it - t') (26) 

V 

the noise becomes thermal, with the strength A = Ao/f , and t)) r]{x' , t')) = 
2T5{x — x')5{t — t'). The QEW becomes then the normal Edwards and Wilkinson 
one, for a surface relaxing by surface tension flOl . 

The QEW develops critical correlations in the vicinity of the critical point, 
Fc- The standard analysis of the problem is the functional renormalization group 
method. One-loop expressions for the exponents are found in papers by Natter- 
mann et al. and Narayan and Fisher l60l l59l . The analysis has been pushed re- 
cently further by Le Doussal, Wiese and collaborators |83| illuminating several 
open issues. The main fixed point is the RF one, but for correlated noise there is 
the possibility of continuously varying exponents. The upper critical dimension 
is dc = 4, above which mean-field theory applies. Several exponent relations ex- 
ist close to the critical point, like 9 = ^{z — x) for the velocity exponent, and 
X + '\-/v = 2. This, together with the ^-exponent relation, tells that there is only 
one temporal and one spatial scale at the critical point. 

For the RF fixed point, the 1-loop functional RG results are x = (4 — (i)/3, 
and z = 2 — {A — d)/^. Later work by Chauve, Le Doussal, and Wiese |83 1 yields 
(with e = 4 — d) C = §(1 + 0.14331e + ...). In particular in d = 1, this means 
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that the RG beyond one loop is better able to adhere to the numerical LEVI results 
EHJI^IHll with a superrough interface (XQEW,d=i ~ 1-2 .. . 1.25). Note that the 
depinning problem can also be discussed in the constant velocity ensemble, which 
allows for the presence of avalanches (since parts of the interface are pinned close 
to Fc) in the presence of translational invariance 

There is one "easy" case in which one can understand SOC with depinning, "at 
once". This comes in the form of the Id Oslo ricepile model (WT\, which operates 
on the idea of having slopes Zi = Ui — n^+i, at each site, which are compared with 
a critical slope Zc, either one or two. The dynamics is such that a overcritical site 
gives one grain to its right neighbor, and thus decreases by 2, while its nearest 
neighbors' increase one. After each toppling Zc is redrawn from a probability dis- 
tribution. If there are no active sites, the avalanche stops and the model is driven 
at the boundary (i = 1). Paczuski and Boettcher noticed that this can be mapped, 
(since Zi changes as an elastic force term) into a boundary-driven LIM |5|, and 
conjectured that the discrete equation is in the LIM class. Recently Pruessner has 
shown that one can do the mapping in a slightly different way, which allows to get 
around the delta-functions involved in this particular model systematically L88J . 

4.2 Mapping SOC to interfaces: the noise 

One would like to describe sandpiles with the dynamics of the "interface", of the 
history of the sandpile, as in the discrete version of a QEW of Section II, the 
Leschhom automaton. For this purpose it is also instructive consider some other 
models than the "LIM" sandpile, the BTW and the Manna ones [81 . 

The Zhang model fT6l resembles both the BTW model and the Manna model 
in that the topping rule is a combination of deterministic and stochastic factors. It 
is defined in with a continuum 'energy' z, with the dynamics 

z{x,t + l) = 0, (27) 
z{y,t + l) = z{y,t) + z{x,t)/2d, (28) 

and Zc = 2d — 1. These imply that the 'energy' of a critical site is divided equally 
between the neighbors. The drive in the Zhang model is usually, in the SOC case, 
implemented so that the energy is increased by small, finite amounts at a random 
site X. Or, one picks the site closest to Zc, since it will typically sooner or later 
be the one to reach criticality first. In the Zhang model the detailed history of 
topplings is of importance making serial and parallel dynamics for active sites to 
differ. Note that the Leschhorn (QEW) automaton and the BTW one are Abelian 
(as is the Manna model, in the general sense), so that it explicitely does not matter 
if parallel or serial topphngs are used. 
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Another widely studied model is the Olami-Feder-Cristiansen (OFC) model 
(§3, usually studied in the following deterministic version. The system is prepared 
with an initial (continuum) energy profile z{x,0) and a threshold Zc is chosen. 
Now, define a distribution parameter a. With the aid of a the OFC model reads if 

z{x,t + l) = z{x,t)-Zc, (29) 
z{y,t + l) = z{y,t)+azc, (30) 

with a < 1/4 in 2d . The point of dissipation (a < 1/4) is that one can achieve 
a steady-state even if the open boundary conditions are replaced with e.g. periodic 
ones. The same trick can also be applied to any of the previous models, where one 
can do away with dissipative boundaries, and substitute for that with a small but 
finite removal rate of grains that topple (HI- 

A third case, very similar to the Manna, is a microscopic model in the NDCF 
class |63|. Consider a two-species reaction-diffusion process, with particles of 
types A and B involved. At each site i, and at each (discrete) time step the follow- 
ing reactions take place: 

Bi ^ Bnn, rdi= 1) (31) 
Ai + Bi ^ 2xB,, n (32) 
Bi ^ A, r2. (33) 

The Ai, Bi denote particles of each kind at site i. r's are the probabilities for the 
microscopic processes to occur: diffusion, r^, activation ri, and passivation r2. 
Fix Vd = I, implying that, after having the chance to react, B particles diffuse with 
probability one, and a phase boundary follows between the active and absorbing 
phases in terms of the ri, r2 probabilities. For occupation numbers nA,i, nB,i per 
site, since the A's are non-diffusive, there is an infinite amount of absorbing states 
defined by ns^i = for all i, with j arbitrary, similarly to the Manna model 

Ell- 

The idea of how to map the dynamics of sand in these models into a noise 
term is similar to the microscopic arrival-time mapping which connects directed 
polymers in a random medium to the temporal behavior of a roughening interface, 
governed by the Kardar-Parisi-Zhang equation |10|. It is also a cousin of 'run- 
time statistics' by Marsili et al., which maps the quenched disorder in say invasion 
percolation to an effective memory term at each location ll90l . The ensuing noise 
in a sandpile LIM arises from annealed to quenched disorder mappings. 

Essentially, one constructs the "local force" (/) by looking at how z^ deviates 
from its expected value at the time of toppling. This is obtained by considering the 
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toppling 



Figure 3: How the randomness maps to an interface noise field in the Manna model 
in Id: each toppling implies one expected grain at a nn, but can give zero or two 
instead (contributing to 6n^^n)- 

net effect of the sandpile rules, by decomposing via a projection trick the expected 
incoming grain or energy flux as 

= njr + <5njr, (34) 

where n^^ measures the expected flux into site x up to time t and 6n^^ is the devi- 
ation. The average part can be used to construct a Laplacian, while (5n™ generates 
a noise term 

t{x,H) = 5nf = nf-nf 

= ni^-uY,H{xnn,t). (35) 

The sum over the nn's in Eq. is the average grain flux into x due to each 
toppling of a nn, computed exactly at the time of toppling at x at height H. For 
the reaction-diffusion model one needs to define H by the integrated activity of 
S-particles, so that one B diffusing out of x means a "toppling". Then the flux 
fluctuates since B's diffuse, randomly, while the reactions between the A and B 
species are accounted for every time B particles leave the site x. le., one con- 
siders riA and ub when the site becomes active and a particle diffuses out. Then 
ntotix, H) = UAix, H) + nB(x, H). 

In the case of the Manna model, each nn-toppling contributes 1/d to and 
thus u = 1/d while for the others the average flux is dependent on the average 
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energy at the time of toppling, 

6n}^ = n}^-Za.,cY.Hnn. (36) 

nn 

Here Zav,c is the above average energy in the Zhang and OFC models, contrary to 
the Manna model this is not known a priori. Another way to write the fluctuation 
term is (5n*" = y^,u Szm , ie. the fluctuations consist of sum of the differences 
between the real amount of energy of the toppling neighbors and the expected 
average energy. This makes it evident that though the Zhang and OFC model will 
have the same columnar noise from the drive as the others, another component is 
born out of the the sandpile rules. 

In the RD-model, the noise takes into account that since after a toppling a 
passive site can get activated only after a nn topples, and since the reactions make 
ub fluctuate while it is non-zero. Thus if and only if a. site does not topple during a 
time step, it implies that there is an effective threshold so that ub — ric < or that 
Ub = 0. These give rise to a r-distribution, that depends upon the total number of 
particles after the preceding toppling and the microscopic dynamical rules at a site. 
The immobile grains ua imply a "pinning force", as a large ua for a constant ntot 
reduces the probability to topple. 

The noise r(x, H) is dependent on the choice of dynamics or the toppling 
order, and this becomes important if the model is not Abelian (the Zhang model) 
Keep in mind, that the mapping describes any particular sandpile with an interface, 
which follows exactly the same dynamics. For any change of dynamics this implies 
that the noise field has to change as well. The / of the discrete interface equation 
reads thus 

f = u V^H + F{x, t) - Zc{x, H) + r(x, H) (37) 

allowing for varying Zc for the sake of generality, with v = \/dfor the Manna. 

The other source of "SOC" noise arises since the implied step-function, 9{f), 
in Equation ^ forces the corresponding interface to move only forward, so that the 
velocity v = AH/ At is either or 1 |58|. This can also considered as an origin 
of noise, a term a{x,H). An illustration is shown in Fig. |4] On the avalanche 
timescale / < at site x and / increases until at time t — 1 one or more neighbors 
topple resulting in / > 0. Then site x topples at time t. Thus the sandpile rules 
result in an effective force /' = 1, acting on H as AH /At = f 9{f) = f'e{f'). 
The relation between / and /', valid for as long as H is constant, is 

f'{x,H) = f{x,H)+a{x,H), a{x,H) = I + z,{x,H) - z{x,t*) (38) 

where t* is the time at which site x topples. This means specifically that one can 
construct cr(x, H) as a quenched random variable, at t* , as if it were noise included 
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Figure 4: Rescaling of the force / and an example of how the cr-noise ensues (three 
grains added simultaneously) due to "extra" grains present when the site topples. 

in a QEW from the very start. It is computed from the difference between / and /' 
when X topples. /' and / are time-dependent, they change as grains are moved, or 
the Laplacian changes. The equation AH/ At = f'9{f') can finally be viewed as 
the discretization of the continuum equation dH/dt = /'. 

4.3 Sandpiles at cridcality 

Such mappings result in the discretized interface equation 



with the surface tension v depending on the model, and with ri{x, H) = —Zc{x, H)+ 
t{x,H). This is the central difference discretization of the LIM, and parallel dy- 
namics are usually assumed. The physics of such an equations contains the Lapla- 
cian character of grain dynamics, which corresponds to the elasticity, topplings 
map exactly to an elastic force. The sandpile dynamics (a toppling) is translated 
into a force rj + F, that manages to overcome a pinning force, so that the interface 
moves by one step. The rules of the individual sandpiles are embedded in the noise 
variables [F, r]{x, H)] as do the details of the dynamics [a]. The ensemble (SOC, 
FES or depinning etc.) is reflected in the boundary conditions, eg. H = at the 
boundaries. The character of the criticality is therefore determined by a combina- 
tion of these two: noise, and ensemble. 

The theoretical understanding of the relevance of the noise terms implies that 
pertubations of the type r/ are relevant, and lead to the establishment of unversahty 
classes that depend on the RG flow of the noise field upon rescaling. In particular, 
there is the LIM/QEW class. As discussed earlier the NDCF/Manna class defines 
another one, characterized by a different set of exponents. The upper critical di- 
mension of the LIM is (i„ = 4 for all cases, however, because of the Laplacian 



dH 
dt 



(uV^H + r]{x, H) + F{x, t) + a{x, H)) 
X e (vV'^H + ri{x, H) + F{x, t) + a{x, H)) 



(39) 
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in Eq. ( l39l . Table 2 lists some recent numerical results of Chate and KockeLkoren 
1^ . The BTW model presents a different story, due to the fact that the continuum 
and discrete versions differ (as we discuss below in the context of the cr-term). 
Notice that its continuum version is directly solvable, and corresponds to the deter- 
ministic relaxation of an initial height profile in the absence of noise, of any type. 
This is since the height profile can be recast into a form that accounts for F{x, 0). 
For the OFC and Zhang models the situation is not so clear - they both present de- 
terministic dynamics during the avalanches, ie. an absence of explicit noise terms 
of the type Zc, but F{x, t) is still present. 

While ri{x,H) term in the LIM case has an exact equivalence Zc{x,H) = 
r]{x, H), implying point disorder, the noise field for the NDCF class, including the 
Manna model has a noise field t{x,H) which is point-like but correlated. This is 
since r is conserved, 

J2 t{x,H) = (40) 

x,H{x) 

in the case of a closed system (periodic boundary conditions, e.g. the fixed energy 
ensemble), and approximately for the SOC one. The average of r is zero, (r) = 0. 
For the Manna model at any particular x the increments of r are random walk -like, 
since the neigboring sites' topplings constitute a random point process. Define the 
noise two-point correlation function by 

C{Ax, AH) = {t{x + Ax,H + 5H)t{x, H)) , (41) 

where an average is performed over the 'disorder' , or in the sandpile sense many 
toppling histories. The random walk-character of the r means that 

C{Ax, AH) ~ f{Ax){AH)^/^ (42) 

since the increments are uncorrelated at a singel x, and / denotes a correlator in 
the x-direction. This kind of noise correlations are similar to the ABBM-model 
of Barkhausen noise. The character of the x-part of the noise correlator for the 
Manna-model is more subtle (note that one might also study such a LIM, but with 
t{x, H) changing after each toppling with uncorrelated increments). In the SOC 
ensemble the interface will be parabolic, in such a way that at a constant H the 
various noise values t{x,H) decorrelate so that one should perhaps measure the 
noise correlator along the interface. Since it is the sandpile dynamics that gives 
rise to r, its local value should be related in a non-trivial way to the local interface 
height. Even more simply, if a site gets by chance very few grains compared to the 
expectation, then it does not topple very much {H < H) and t(x, H) is negative. 
In the FES case / can be defined and measured e.g. numerically in the standard 
fashion, but C ^ f should decay quickly since neighbors have a weak influence on 
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each others' noise values r (see Fig.|3ll. They correlate in the noise strength since 
the H-pait of the correlator is dependent on the interface height. The case of the 
RD model is slightly more complicated, though the noise from the diffusive flux is 
the same as for the Manna, explaining qualitatively the existence of such a NDCF 
universality class. The off-critical behavior of such models couples directly the 
noise field, unlike in a LIM-automaton, and might have to do with the observation 
of slow relaxation in the Manna model |92 1. 

Keep in mind that what matters in determining avalanche properties is the ef- 
fective disorder that influences its spreading. In the models with r-noise this means 
that increments of r stem from randomness in the rules (and the noise correlator C 
above measures its feedback on the avalanche dynamics). For other models (BTW, 
Zhang, OFC without additional disorder) the avalanche dynamics does not experi- 
ence any randomness "during" the avalanche, but are set by the initial properties of 
the system. Eg. though for Manna r is conserved, in the case of the Zhang model 
(say) the average r during an avalanche can be different from zero, though its en- 
semble average is zero. It is also possible to drive a SOC-interface with strong 
enough randomness so that the interface tension vanishes. In sandpile language 
this can be achieved by making e.g. the toppling thresholds increase (even without 
bounds). Then a natural possibility is a cross-over to directed percolation, precisely 
since the diffusive grain dynamics is changed, by the extra noise |93|. 

The fixed energy ensemble |13, 68 1 corresponds to ordinary depinning, ie. 
H{x,t = 0) = 0, F{x,t) = F{x,0), and periodic BCs, with no dissipation. 
One needs to tune the control parameter AF = ritot/L'^, to obtain criticality, to 
A = 0. The usual "random preparation" of the grain configuration corresponds to 
a spatially dependent force profile F{x, 0), which might give rise to some memory 
effects, depending on its exact (columnar) form l94ll . In the same vein, "micro- 
canonical" simulations [951 one has dissipation operating on the slow time scale 
with exactly the same rate as F{x,t). Thus microcanonical simulations corre- 
spond to fixed energy simulations with a specific initial configuration: after each 
avalanche, the time is reset to zero, the force is replaced with F —>■ F -\- V^H, 
and the forces at x' (x") are increased (decreased) by one unit where x' and x" 
are randomly chosen sites. Finally the interface is initialized, H = 0. Recent 
simulations of the NDCF/Manna (or Conserved-Directed Percolation, C-DP) class 
in the depinning ensemble have revealed, that the scaling exponents are in 2d 
very close to LIM, while in Id and 3d differences ensue [91 96, 65 1. In partic- 
ular, there is the K-exponent that relates the local and global roughness exponents 
(Xgiobai = 2k + xioc, as a clear signal of anomalous scaling (and in Id of super- 
rough behvior) 19111961 : the . The BTW model is in a different class as such, due to 
the extra symmetries, and as expected based on the presence of noise in the other 
models |97|. 
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Table 2: Critical exponents of the C-DP/NDCF class in Id, 2d, and 3d. The 
corresponding values for DP and LIM are also given for reference. From L91J . 



d 




6 


z 


P 


2k 




DP 


0.1596 


1.58 


0.2765 


0.84(1) 


1 


C-DP 


0.140(5) 


1.55(3) 


0.29(2) 


0.86(1) 




LIM 


0.125(5) 


1.43(1) 


0.25(2) 


0.35(1) 




DP 


0.451 


1.76 


0.584 


0.56(2) 


2 


C-DP 


0.51(1) 


1.55(3) 


0.64(2) 


0.50(2) 




LIM 


0.50(1) 


1.55(2) 


0.63(2) 


0+ 




DP 


0.73 


1.90 


0.81 


0.30(5) 


3 


C-DP 


0.88(2) 


1.73(5) 


0.88(2) 


< 0.2 




LIM 


0.77(2) 


1.78(7) 


0.85(2) 






For open boundary sandpiles, the right interface boundary condition is H = 
which is to be imposed at "extra sites" (x = 0, x = L + 1 for a system of 
size L in Id). In the SOC steady-state the Laplacian increases, because of the 
parabolic shape for the the toppling profile H{x) (see also l98l : p in the APT- 
picture naturally equals the same). This is compensated by the ever-increasing 
{F{t)), or the addition of grains by the SOC drive. Writing now (F) = ft, and 
H = Hxt, one has the Poisson equation 

uV^Hx = -U (43) 

with the constraint that H = 0, x € dx, and the source term fx = 0,dx and 
fx = /> elsewhere. This formulation of the 'history' of topplings in a sandpile is 
just another way of looking at the lattice Green's function for any sandpile, noticed 
by Dhar in the case of the BTW long ago |99|. 

Concerning the fluctuations, at the depinning critical point an interface has a 
diverging response. In the sandpile language the average avalanche size diverges 
with system size, 

(s) ~ (44) 

which follows also from the fact that each added grain will perform (an effective 
random waUc) of the order of topplings before leaving the system, independent 
of dimension [13 100|. In Eq. (143 1 the rate of divergence is reflected in u, in 
that it measures for each particular choice of rules how the above scaling can be 
interpreted as an equality, since the total number of topplings per grain added is 
dependent on v. The linear relation between v and the interface velocity does in 
no way reflect the actual universality classes of the models. It just manifests the 
fact that the critical point of such models exhibits "anomalous diffusion", as hinted 
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by the LIM dynamical exponent z ^ 2. Hence in the interface representation such 
glassy response is easy to understand IIIOIII . 

Recall that in the SOC steady-state the probability to have avalanches of life- 
time t may have the distribution P{t) = t^'^^ ft{t/L^). P{s) was mentioned in the 
Introduction, and s ~ t^^/^ and z{Tt — l) = (r^ — 1). The linear size £ related to 
s scales as s ~ £^ and the (spatial) area as i'^ (for compact avalanches) with £ the 
linear dimension. Thus, t = 2 — 2/ D and Tf = 1 + {D — 2)/ z. From the interface 
depinning picture we now obtain directly that the "volume" of an avalanche should 
scale as d+x^ ^^'^ that the average area scales as {^^"^^ ~ L^/*^. The former implies 
the important scaling relation 

Ts = 2-2/D = 2-2/{d + x)- (45) 

This holds if and only j/ simple self-affine geometry can be expected. The avalanche 
exponents of SOC models are notoriously hard to measure, making it for a long 
time apparent that e.g. the BTW and Manna models would be in 2d the same 
universality class, with ~ 1.27 for instance llT02l fT03l fT04l . In the light of 
connection to interface depinning and (Manna, only) to NDCF class this seems 
fortuitious at once. The constant-velocity ensemble is a "good" one for checking 
out Eq. (l45l . and indeed for the LIM class the expected avalanche scaling can be 
found 1 86 1). 

For the SOC case, recent simulations have unearthed two fundamental facts. 
First, the dissipating avalanches can have interesting properties of their own 1 105 1. 
In the BTW model this is related to the fact that the avalanche stops at once when 
the boundary is met, while for the other models this holds as a strong trend i22l. 
More importantly, the existence of multifractal behavior has been demonstrated, 
also in the BTW, by Stella and coworkers f24|. This presumably follows from the 
strong symmetries present in it; one can also measure many interesting quantities 
in the "wave picture" l2Tlll06lll07l . of the BTW, which are naturally hard to relate 
to any scaling theory related to interfaces or absorbing states. In concrete terms, 
the multiscaling indicates that for e.g. the avalanche sizes various momenta exhibit 
different scaling exponents Tg^q, as a function of the pile size L, {s'^) ~ Ll ^. For 
the Zhang model, Pastor-Satorras and Vespignani have used the same technique 
to point out the lack of any clear scaling whatsoever. This model can combined 
with randomness, in which case the scaling attained is that of the NDCF/Manna- 
class. Likewise it seems possible to pertub a SOC state continously (in terms of the 
effective scaling exponents) between the Manna and the Zhang model endpoints 
II108I . This seems to follow naturally in the light of change in the effective noise 
acting on the interface representing the sandpile. In the same vein, the existence 
of strong. Manna-type noise is a relevant perturbation and explains the cross-over 
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to similar scaling exponents, if such randomness is added on the top of the other 
sandpile rules 1 109 1. 

Recent work by Dickman has shown however that more fundamental surprises 
can be found in the NDCF avalanche properties 1 1 10 1. Fig. Q depicts the result of 
careful numerical simulations, demonstrating that one needs systematic logarith- 
mic corrections to describe the P(s)-distribution lllll . That is, 

P{s) = s-^' {In syfs{s/sc), (46) 

where r = 1.386 7 = 0.683, for the Manna model in 2d (Id data did not show 
evidence of such corrections). /* (see the caption) fluctuates about a constant over 
the optimum fitting interval, while the effective r varies strongly in the same. The 
figure also has a pure power-law fit, with the estimate = 1.25 11031 . The latter 
yields a strongly curved function and shows that a power-law is not a good 

description. 
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Figure 5: Plot of /* = s'^ Ps{s) / {In s)'^ versus Ins for the data shown in Fig. 1. 
Lower curve: best-fit for 3 < Ins < 10 using r = 1.385 and 7 = 0.672; upper 
curve: pure power-law fit using r = 1.25. From lllll . 

The next candidate, to establish the validity of Eq. (l45l . is then the LIM model. 
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It also allows to use the geometry of the "pinning paths", if the SOC and normal 
critical states are directly connected. In directed percolation depinning (DPD), 
the geometry of the random quenched landscape leads to a self-affine picture of 
the interface dynamics, described by the quenched KPZ equation |72, 10 1. The 
interface motion takes place so that the interface invades the voids (constituting 
avalanches) of multiconnected network of pinning paths, and thus the sizes of the 
avalanches are related to the void sizes. This is given by (in the LIM case) the 
structure of an "elastic percolation problem" 1 112J . so that the the RHS of the qEW 
is always negative semi-definite, / < 0. There are two fundamental issues: how 
the scaling of voids, the relation of size vs. area is described with an appropriate 
roughness exponent, xioc, and what is the probability to produce an avalanche if 
the interface is unpinned at a particular spot ( 11131 II 141 III5I ). In analogy with 
DPD L113J .it follows that 

Ts,dep = 1 + (1/(1 + X/oc))(l - l/iy) (47) 

which using the global x ~ 1-20 ... 1.25, produces 

Ts,dep ~ 1.08. (48) 

This is the same as the prediction of Eq. (l45l . showing that assuming strict self- 
affine avalanche scaling the SOC and depinning ensemble predictions for (and 
other such exponents) coincide. 

L= 128 ... 8192 

10° I I I I I I I I I 




^g-6| I I I I I I I I 

10"'° 10° 10° 10"' 10"° 10° 10° 10* 
s/L° 

Figure 6: A data collapse using ~ 1.08 and Ds ~ 2.20 for the LIM-SOC model 
(Leschhorn-type automaton, unpublished). 

Figure|6lshows a representative scaling plot of SOC simulations of a Id qEW/LIM, 
of the Leschhom-type. System sizes upto L = 8192 have been studied, with 2 x 
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10^ avalanches for the largest L. It is rather evident that the bulk part of the dis- 
tribution does not collapse that well, for the largest system size the fitted effective 
T is way off from the predicted 1.08 (1.024 for L = 8192). An attempt to fit the 
exponents to an asymptotic value and a correction implies Ts^soci'^) ~ 1.115, 
clearly off the expected value. 

Thus, in analogy to the 2d Manna model, it seems that the ensemble here is of 
importance, or in the interface language the parabolic interface profile. It would 
be of interest to attempt such comparisons in higher dimensions, where one needs 
more exponents (or assumptions to describe the probability distributions since the 
avalanches have in addition to an area vs. volume relation also a perimeter length 
vs. area relation 11161 11141'). One would expect in any case that < d + x for 
the cut-off dimension, since the SOC state can of course not be 'over-critical'. 

Surprises in trying to match the translation-invariant non-equilibrium phase 
transition and the SOC one present the question whether the differences vanish 
asymptotically, so that the exponents of the statistically homogeneous ensemble 
are recovered? The other possibility is that the SOC ensemble is an independent 
one. Microscopically, this would result from a non-uniform density of grains (or 
average force acting on the interface). In Id, a site x will get a larger grain flux 
from its neighbor on the bulk side and a smaller one from the boundary side - more 
trivially, there is a net flux of grains towards the boundaries. 

This can be stated a bit less casually by considering the integral of the deviation 
from the force at depinning criticality, at x, AF{x) = F^ep— F {x) soc ^ where F^ep 
and F{x)soc are averages at the critical points of the ensembles. This implies e.g. 
a finite-size correction to the normal critical point, / /S.F{x)dx = 5F{L). 6F of 
course follows from the exact scaling function of F{x)soc- In analogy to normal 
depinning, this defines a correlation length exponent i^soc via 5F{L) ~ ]J-/'^soc_ 
Usual critical exponents like v are derived from the RG, assuming an ensemble 
with statistical translational invariance. In the SOC case this is lacking, and it is 
not immediately clear whether one can obtain the properties of the SOC state from 
e.g. boundary criticality II117I . In this respect the usual SOC models are more 
complicated than e.g. the boundary driven Oslo model. 

Fig. Q shows, in contrast, a counter-example, from the Oslo (ie. a boundary- 
driven LIM) model. Here the continuum equation (of course integrated numeri- 
cally) obtained from a mapping of the sandpile to the LIM is compared to the usual 
automaton for the avalanche size. The result is that the exponents coincide, with 
each other and with the analogous prediction for boundary-driven avalanches, to 
Eq. d45l (that uses (s) ~ L). This examplifies the fact that the SOC state can be 
understood, if translational invariance is present. 

Finally, the properties of the SOC state are also visible in the effective noise 
terms of the interface (e.g. r). It is not to be expected a priori, that the noise corre- 
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Figure 7: Data collapses for system sizes between L = 128 and L = 512, for 
the continous and discrete boundary-driven LIMs. The same value of = 1.55 
suffices for both 1 88 1 

lator and strength are translation: consider the cr-noise, which reflects the tendency 
of the interface to move faster. This effect is the strongest in the center of a SOC 
model, while with periodic boundaries the strength of the noise is independent of 
X. The origin of a depends on the exact sandpile rules, in the Manna-model a site 
with Zx = I can also get two grains from the same neighbor, making it overactive 
(note that one can thus simply combine a and r to a single quenched noise term), 
while in the Zhang, OFC and BTW models the cr-field is set from the onset of an 
avalanche, including the toppling order chosen. The probability for a{x, H) ^ 
is related to the fluctuations in the density of active (and critical, = Zc) sites 
and is thus not accessible by graph theoretical 1 106 II or mean-field analysis 1 14 1 
of non-active configurations. It is weaker in d > 2 due to the ramified avalanche 
structure (and has no role with respect to the upper critical dimension). Its strength 
evidently changes off the critical point (since more of the neighbors are likely to 
be active at the same time). 

Simulations of cr-fields in 2d reveal that the average noise strength depends on 
L, in terms of Pl(o" < 0), the probability to have a non-zero value at a site at a 
toppling, and (o"(x)) is non-uniform in x. The asymptotic BTW values are {P) ~ 
0.081 for all and {P)^ — 0.121 for dissipating avalanches. Thus the dissipating 
avalanches have typically stronger cr-noise. In the BTW model the noise field 
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Figure 8: P{a{x) < 0) vs. x scaled with P{x = L/2) for various L in the BTW 
model (L = 32 (circles), 64 (triangles), 128 (squares). Inset: for L = 64 the same 
for BTW, Manna, ricepile (triangles, circles, crosses) |8l. 

has exponentially decaying correlations for any fixed x, but with a decay length 
that is largest in the center of the pile, and increases in general with L. This and 
the existence of periodic oscillations in the correlations in the bulk underlines the 
missing translational invariance citedrossel,ktitarev,barrat. It is unclear how the 
correlations of a relate to the multiscaling in the model 1241 . but certainly they 
reflect also the columnar noise of the interface equation II118I . In contrast, for the 
Manna or LIM automata the correlations in a decay rapidly. Figure [8] depicts a 
data collapse of P{a{x,y) ^ 0) along the cut y = L/2, 1 < x < L/2, scaled 
with P{x ~ L/2), evalued for those topplings that are the last at the site during 
an avalanche. The probability increases at the boundaries with system size L. This 
shows that the lack of spatial translational invariance is a generic property of the 
SOC ensemble. 

4.4 Off-criticality in SOC or interface models 

For interfaces the outcome of driving the system "too hard", F > Fc (or preparing 
a sandpile to be "overcritical", with ntot), is just a cross-over to, eventually, thermal 
noise on lengthscales L ^ ^. There are some typical ways to push SOC models 
off the critical point. The first option is to add a loss mechanism, to the dynamics 
of grains, typically in the FES ensemble. The critical state is then reached only 
for carefully tuned bulk dissipation e '-^ L^^ (e.g. |14|). The interface equation 
develops - using the same projection trick as for the Manna etc. models - a linear 
confinement term 

— n^V^H + i]{x,H) -eH{x) + F{x,t) (49) 

where e > measures the strength of the dissipation. The —eH increases with a 
(small) probability only when a site topples and the fluctuations in the dissipation 
contribute to ri{x, H). For e.g. the BTW model this means a perturbation, which is 
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irrelevant as long as the Larkin length associated with the cross-over from columnar 
behavior is larger than L, and thus the avalanche behavior is governed by the BTW 
dynamics. This is in fact the critical point of the constant velocity ensemble of 
depinning. For all such sandpile models one indication is however clear: for large 
avalanches the critical behavior is cut off by the linear term, and thus any value of e 
is sufficient in the TD-hmit to bring the system off the critical point (regardless of 
the ensemble) fll9l. This is due to the avalanches becoming d— 1-dimensional, or 
that any site topples at most once. This result is of relevance for the long-standing 
discussion about the existence of true criticality with dissipation in the OFC-class 
of models, but does not provide any hint about any remaining scaling for a < 1/4 
flM . 

The driven sandpiles with boundary dissipation have been also analyzed by var- 
ious means. E.g. Kardar and Hwa attempted to describe the grain flux dynamics 
by a RG treatment, to understand the correlations in the ensuing steady-state (for 
instance by its power-spectrum) |19|. The interface picture implies a cross-over 
between quenched and thermal noise, which among others means that the interface 
velocity can be related to the fluctuations in the grain flux. Off the critical point the 
avalanches overlap, and the stoppage of activity becomes a rare event, like the in- 
terface velocity becoming zero when F > Fc [121 1. The expectation would be that 
one obtains 1 //-noise in the dynamics |l l22lll23J . The correlations in the activity 
or interface velocity reflect the Laplacian nature of grain dynamics [98). Figure |9l 
shows an example of the response functions of constantly driven BTW model, to 
an extra perturbation. It exhibits clearly the properties of a diffusive response to the 
extra drive field: the response function (/*~°° 5Hdt) decays exponentially with r, 
distance from the location where the perturbation was applied. 

With the SOC criticality destroyed, the (random) drive and boundary dissipa- 
tion still affect the fluctuations. Usual measurements fT9l l98l IT24 I) concentrate 
on instantaneous quantities like the local force/grain density or activity/velocity. 
With a boundary condition H = imposed, a drive {F{x,t)) = ft with / a 
fixed constant produces a constant average velocity which varies with x, {v{x)) = 
{dtH{x, t)). The continuum equation version is 

= uV'H{x, t) + ft + 6 fix, t) + r,{x, H{t)), (50) 

a depinning ensemble with a constant drive rate f . This is not equivalent exactly 
to the normal ensembles (constant force, or constant velocity [85J ). The fluctuat- 
ing part of H, 5H, should reflect the fact that the noise rj{x, v{x) + 5H) develops 
temporal correlations that depend explicitely on x. Translational invariance is thus 
absent also off the depinning critical point. Note that the boundary regions are 
closer to depinning, and lim^^o, L ^^(2;) = 0^, while {v) oc f jl?'. In the prox- 
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Figure 9: Time integrated response function Xh-*Q,e to a constant perturbation as a 
function of r, the strength of the perturbation, the linear lattice size is L = 200. 
The lines are exponential fits. From ref. I.98J 

imity of the SOC critical point one feature of Eq. (l50l seems to adhere to normal 
depinning: the waiting times {v{x) = 0) of the interface follow roughly a power- 
law distribution P{tw) ~ i^^^^™ II125I . where (3^ is the roughening exponent of 
the FES/depinning problem, with no signs of logarithmic corrections 11111 in the 
case of the Manna model 1 126 1. 

The 6H for a Id interface is an example of retum-to-zero stochastic processes 
recently analyzed by Baldassari et al. II127I . Here it is complicated by the struc- 
ture of rj{x,t). For the truly thermal EW equation, without noise correlations the 
average fluctuations form a parabola (in amplitude), which is easy to understand 
since the EW interface profile is a random walk with the given boundary condi- 
tions. Measurements of the amplitude {6'^H)x as a function of x for a Id LIM 
SOC model, driven off the critical point, imply that the function changes from the 
EW form. Likewise, the two-point temporal correlation function of the interface 
has effective scaling behavior that combines the effect of "kicks" (drive by grains) 
and the relaxation of the interface, and also reflects at the boundaries the fluctua- 
tions in the grain flux, from the interior of the system. All such details are awaiting 
real analytical understanding. The fluctuation profiles in Fig. [TOlshow the same 
scaling behavior for all the system sizes as a function of x for the Id QEW: for x 
small {d'^H){x) ~ x^-^, instead of ~ x^ as for a normal EW. 
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5 SOC in the Kardar-Parisi-Zhang depinning 



The second fundamental equation for depinning problems is the quenched Kardar- 
Parisi-Zhang one 1 129 1301, 

5fF(x, t) = vV'^H{^, t) + \{\IH{^, t)f + F + T]{x, H{x, t)), (51) 

A measures the strength of the celebrated KPZ nonlinearity, proportional to the 
squared local interface slope 11311 . The previous analysis of reaching a non- 
equilibrium critical steady-state can now be repeated in the presence of the A-term, 
to study the difference between the QKPZ and the EW equations (a non-linear 
Langevin equation and a linear diffusion equation) II132I . The dynamics of the 
latter is insensitive to initial profiles (unless no other noise is present), while the 
QKPZ has a growth component dtH perpendicular to the interface slope, always 
non-zero. 

The Eq. dSTT i suffers from the problem that it is notoriously hard to integrate 
numerically. Some discretizations exist, a reasonable example 11331 is given by 

m = uV^H^t) + ^ {[Hi{t) - Hi^i{t)f + [Hi+i{t) - H,{t)f) 

+F{t)+r],^Hr (52) 

where fi{t) is again a local force. The QKPZ sandpile is such that a site gives 
2v particles or units of energy to its neighbors once it topples, and the toppling 
criterion also compares the integrated activity at site i to its neighbors. In the 
interface picture Hi{t + 1) = Hi{t) + Q{fi{t)). With this definition, H{x,t) = 
/* p{x, T)dT maps the discrete QKPZ equation to the 'sandpile' and vice versa. 
Now the simple dynamics dtH = 6{f ) and pinning boundary conditions H = 
at the boundaries imply a SOC ensemble, if F{t) is ramped slowly as in usual 
sandpiles. Note that the exponents of the depinning transition depend on e.g. if A 
approaches zero with AF = F — Fc —^ (One can map a directed sandpile to the 
KPZ equation, without considering the history fl341). 

The roughening, the development of average local slopes, depends on the sign 
of A. If F{t) > the interface becomes unstable if A > 0. If A < one can study 
a separation of time-scales: the system is driven so that F{t) is increased after an 
avalanche has relaxed. Fig.^Jshows a series of snapshots from a simulation of the 
QKPZ equation with 7?j i^. = ±g with equal probability. This describes random 
toppling thresholds that change after each toppling in the sandpile language. 

The shape of the interface is triangular in 1-i-l-D as expected with A < (see 
also [135]). For a fiat slope the surface tension is negligible, and a slope arises 
from the balance of the nonlinearity and the driving force, (|V/i|) = \/2C,t/\. 
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This is analogous to the development of the steady-state grain configuration in 
usual SOC. In the ensuing phase usual avalanche properties can be measured, eg. 
Ts = 2.5 ± 0.1, implying that the average avalanche size is independent of L, 
contrary to usual sandpiles with anomalous diffusion, since > 2. 

The avalanches are in addition time-translation invariant only on the average: 
the largest avalanches happen when the roughness of the interface is the largest. 
Local slopes, i. e. \Hi — Hi^i \ alternate between n and n + 1 where n is a positive 
integer. In spite of such deterministic trends there is no characteristic avalanche 
size. Such dynamics (for the size, duration, and support) of avalanches eventu- 
ally changes due to the development of a critical slope rric. This is a signature 
of different effective dynamics of avalanches, above nic, which one reaches in 
avalanches. 

The combination of slow drive and boundary pinning eventually pick a sym- 
metry for the average interface profile. Write the interface field on one of the sides 
of the triangle (Id) as H{x, t) = m{t)x + 6H{x, t). Inserting this into the QKPZ 
equation we obtain for 6H 

^'^^i'"'^^ = uV^6H + Xm{t)V6H + ^{V6Hf + 7?(x, H), (53) 
at 2 

subtracting the mean-field solution, valid for intermediate times 6t <^ t. The ef- 
fective interface equation has a linear term in V6H, that will dominate over the 
KPZ-nonlinearity. This is a full analogy of depinning in the presence of anisotropic 
quenched noise 11361 . explaining the presence of this regime. 



6 Conclusions and discussion 

In this article we have outlined the connection of sandpile-like cellular automata, 
exhibiting "self-organized criticality", to other non-equilibrium phase transitions. 
Such mappings allow to establish a number of facts about sandpiles, and leave a 
number of open questions. 

The first observations concern the role of ensembles and rules. The mappings to 
interface depinning reveal, how SOC avalanche properties depend on the details of 
the SOC automaton. The connections to both absorbing state phase transitions and 
the depinning illuminate the role of history effects in sandpile dynamics. Below, a 
few evident conclusions are listed. 

1 . the noise in the LIM-representation settles the universality class, in the ab- 
sence of strong additional symmetries (BTW model). 

2. the upper critical dimension, in the presence of a surface tension, is four. 
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3. for non-SOC ensembles the critical exponent follow from the FES/depinning 
transition. In some cases these present an open problem, since the QEW 
with certain kinds of correlated noise is not well understood (in terms of x, z 
etc.). This has not been resolved in the NDCF-picture either, due to the lack 
of analytical progress. 

4. the derivation of the avalanche exponents for SOC sandpiles is an open prob- 
lem. In some cases (Zhang, Manna/NDCF class) this can be directly traced 
to the fact that the noise of the interface equation becomes x-dependent. 
Thus translational invariance is broken. The same holds for off-critical SOC 
sandpiles. 

5. the off-critical models can still show criticality (eg. of the thermal EW-class) 

6. temperature can be introduced easily II137I . via the idea of adding "thermal 
noise" into the interface equation II138II1391 . and interpreting it in the sand- 
pile language. 

7. the rule of conservation laws is highlighted (Laplacian operator, the FES 
class and its description via an activity picture). 

8. the boundary conditions of the SOC ensemble may (QKPZ "sandpile") lead 
to different dynamics, due to symmetry-breaking. 

In summary, the existence of "SOC" is by such mappings revealed to be a result of 
combination of "non-linearity" (underlying ordinary non-equilibrium phase tran- 
sition) together with a slow drive in an inhomogeneous system (net current of par- 
ticles or energy). It is not an "attractive fixed point", but arises from fine-tuning the 
drive rate so that avalanches separate. 

One outcome is naturally the prospect of observing "true SOC" by preparing 
an experiment with boundary conditions that exactly correspond to what the theory 
would imply. Driven interfaces in random media (domain walls in random mag- 
nets, for instance) would seem to be an obvious candidate. Note that the predictions 
for open systems differ depending whether the system is exactly in the SOC state 
or not. Concerning theoretical aspects, the mapping of dynamics to a history de- 
scription might be of advantage in other problems, like coupled map-lattices (for 
the contact process, e.g. a new exponent arises 1140111411 ') and reaction-diffusion 
systems (eg. 1 142 1). 

For SOC-like cellular automata and their continuum descriptions one may list a 
number of open topics, including those related to the list above. Understanding the 
nature of the SOC ensemble; the origins of multiscaling therein; mapping depin- 
ning to absorbing state phase transitions 1 7 1 1 renormalization techniques for both 
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depinning (in non-standard ensembles like the SOC, and with correlated noise) and 
the conserved field theory of NDCF models. 

Other topics are the continuum descriptions of models with more complicated 
rules (curvature -dependent thresholds, like for sandpiles derived out of driven 
MHD equations 11431 . in both SOC and FES ensembles; the vortex flow model 
of Paczuski and Bassler that may be in the randomly boundary -driven colum- 
nar noise LIM-class f43l, the non-local rules employed already by Kadanoff et al. 
(Too I and so forth. In the interface context, the presence of long-range interactions 
substituting the Laplacian surface tension in Eq. d39l changes the scaling expo- 
nents 11441 . For all such cellular automata one expects that they can be mapped to 
known continuum equations with quenched randomness, possibly with noise cor- 
relations. This also implies that e.g the predictability |145| of such SOC systems 
is very low and that the damage tolerance is high 1,146J. An analogy is provided 
by a percolation problem, in which one removes bonds around pc randomly till the 
spanning cluster is broken; then new ones are added again till a spanning cluster 
is re-formed. Clearly it is to first order difficult to "predict" whether the largest 
cluster is system-spanning or not and in fact the only information hidden in the 
behavior is the usual spanning probablity (vs. p). 
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Figure 10: Amplitude of the interface fluctuations (normalized with SH{x = L/2) 
{6'^H){x) in Id for both the QEW and a thermal EW model with similar bound- 
aries. The drive rate, for grains, / = 1/100 for the former, for all the L = 
64 . . . 256. flM 



900 




Figure 11: Typical interface configurations for a Id QKPZ, with SOC drive and 
boundary conditions, for a system of 500 sites. The one with the flat top is an 
unstable one from the initial growth event before steady-state of the model, and the 
others are from every 20000th time step I132II . 
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